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ABSTRACT 

Starting with the well-known NFW dark matter halo distribution, we con- 
struct a simple polytropic model for the intracluster medium which is in good 
£T) '. agreement with high resolution numerical hydrodynamical simulations, apply 

this model to a very large scale concordance dark matter simulation, and corn- 
ID ! pare the resulting global properties with recent observations of X-ray clusters, 
O • 

6 ■ 
a: 

good approximation to the internal gas structure of massive clusters (except in 
j> | the very central regions where cooling becomes important), and allows one to 

recover the observed M 50 o — T, L x — T and T/nl^ 3 oc T a65 relations. Using 
these concepts and generalizing this method so that it can be applied to fully 
three-dimensional N-body simulations, one can predict the global X-ray and SZE 
trends for any specified cosmological model. We find a good fit to observations 
when assuming that twelve percent of the initial baryonic mass condenses into 
stars, the fraction of rest mass of this condensed component transferred back to 
the remaining gas (feedback) is 3.9 x 10~ 5 , and the fraction of total pressure from 
a nonthermal component is near ten percent. 

Subject headings: cosmology: theory — galaxies:clusters:general — intergalactic 
medium — X-rays:galaxies:clusters 



including the mass-temperature and luminosity-temperature relations. We make 
allowances for a non-negligible surface pressure, removal of low entropy (short 
cooling time) gas, energy injection due to feedback, and for a relativistic (non- 
thermal) pressure component. A polytropic index n = 5 (T = 1.2) provides a 
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1. Introduction 

Gas in clusters of galaxies can be observed to large cosmological distances by a variety 
of techniques, from X-rays (Bremsstrahlung) to radio (S-Z effect). But to utilize these 
observations it is necessary to have a model for the state of the gas in clusters that is (a) 
motivated by sound physical reasons, (b) able to fit the observational data and (c) is simple 
enough (i.e., not a numerical simulation!) to be applied broadly. In this paper we attempt 
to present such a model for the gaseous component in clusters of galaxies in order to provide 
predictions for global properties such as temperature and X-ray luminosity. 

For the dark matter (DM), the widely utilized Navarro, Frenk & White (1997), or NFW, 
model satisfies all of the above criteria. Although we now know that it is not universal in 
two senses — large variance in the central concentrations (Avila-Reese et al. 1999; Jing 2000; 
Bullock et al. 2001; Klypin et al. 2001; Fukushige, Kawai & Makino 2004; Tasitsiomi, et al. 
2004) and trends in properties with time and halo mass (Wechsler et al. 2002; Ricotti 2003; 
Zhao et al. 2003; Weller, Ostriker & Bode 2005; Salvador-Sole, Manrique & Solanes 2005)— 
it remains an extremely useful first basis for analyzing and summarizing the properties of 
dark matter distributions. We know, however that the gas in clusters does not follow the 
dark matter profile. The well established central density profiles for dark matter halos are 
roughly power laws: density p depends on radius r as p oc r~ a , with a typically 1.0 (the 
NFW value) but ranging from 0.5 to 1.5, depending on circumstances. But the gas profiles 
show a definite core (a — > 0; for a recent review see Voit 2004) and, as we shall see, one 
would overestimate the X-ray luminosity by a large factor if one were to use the NFW or a 
steeper profile. 

The construction of a satisfactory model will be guided by a few observed properties. 
First, the gas is essentially a trace component (approximately 1/7 of total mass) and resides, 
in close to hydrostatic equilibrium, in a potential which is well represented by NFW (or its 
variants- cf. Zhao 1996). Furthermore, we know that there are efficient means of redistribut- 
ing energy/entropy within the cluster gas via, for example, turbulence (Kim & Narayan 
2003b) induced by merger shocks (Bryan & Norman 1998) and galaxy wakes (Stevens, Acre- 
man & Ponman 1999; Sakelliou 2000). Additionally, other processes such as conduction 
(Kim & Narayan 2003a; Dolag et al. 2004), cosmic ray transport, and magneto-sonic wave 
transport (Cen 2005) may also be operating. Also, appropriate boundary conditions are re- 
quired, since in both simple analytic models (eg. Bertschinger 1985) and detailed numerical 
simulations (eg. Bryan & Norman 1998; Frenk et al. 1999) the hydrostatic portion of the 
cluster gas is terminated at an outer shock where the pressure is balanced by the momentum 
flux of the infalling gas. 

A number of steps have already been made toward constructing such a model. Makino, 
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Sasaki k Suto (1998) derived an analytic expression for a gas distribution in hydrostatic 
equilibrium with an NFW potential, assuming isothermality. A more general expression for 
a polytropic equation of state was derived by Suto et al. (1998); a similar functional form 
has been compared in detail with hydrodynamic simulations by Ascasibar et al. (2003). 

In the simplest of such models the source of the gas heating is gravitational, i.e. the 
gas energy comes from the same collapse and virialization processes which determine the 
dark matter profile; thus the energy per unit mass in the gas should be approximately the 
same as in the dark matter. This leads to an expected self-similar scalings of mass M and 
luminosity L with temperature T of M oc T 3 / 2 and L oc T 2 (Kaiser 1986). However, this 
expectation is in contradiction with the observed relation, which is steeper (Edge k Stewart 
1991; Markevitch 1998). Kaiser (1991) proposed that non-gravitational energy injection 
could lead to the observed relations. This possibility has been explored in the type of 
analytic model described here using an NFW profile for the dark matter (Suto et al. 1998; 
Wu, Fabian k Nulsen 2000; Shimizu et al. 2004; Lapi, Cavaliere k Menci 2005; Afshordi, 
Lin k Sanderson 2005; Solanes et al. 2005), and with other profiles (Balogh, Babul k Patton 
1999; Babul et al. 2002). The breaking of self-similarity can also be cast as the modification 
of the initial gas entropy by thermal and nonthermal processes, as explored in NFW-like 
potentials by Tozzi k Norman (2001), Komatsu k Seljak (2001), Voit et al. (2002), and 
Dos Santos k Dore (2002); on the other hand Roychowdhury k Nath (2003) argue that 
the entropy imparted to the gas from gravitational processes alone is larger than previously 
thought. Another impact on the gas energy comes from the fact that approximately one 
tenth of the baryons in a typical cluster are now in stellar form. So one must allow for both 
removal of the mass of this gas and of the associated energy (or entropy) of this gas (Voit k 
Bryan 2001; Tozzi k Norman 2001; Voit et al. 2002; Scannapieco k Oh 2004; Bryan k Voit 
2005). Since the removed gas had short cooling times, low entropy, and low total energy, 
the mean energy per unit mass of the remaining gas is higher than before star and galaxy 
formation. 

An issue not dealt with in these studies is non-thermal sources of pressure. Turbulence 
may provide in excess of 10% of the total pressure in Coma (Schuecker et al. 2004); sim- 
ilar amounts of turbulent support have been seen in simulations (Norman k Bryan 1999; 
Faltenbacher et al. 2004). Clusters should also contain a population of relativistic parti- 
cles arising from shocks, as recently reviewed by Miniati (2004), Sarazin (2004), and Bykov 
(2005). Magnetic fields may also be dynamically important (Carilli k Taylor 2002; Ensslin, 
Vogt k Pfrommer 2005; Bykov 2005). 

The basic goal of this paper is to start with a population of dark matter halos from an 
N-body simulation, for which the DM density profiles can easily be measured, and deduce the 
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global properties of the hot baryonic component in a physically well-motivated manner. The 
ideal method should be as simple as possible while including all the relevant components: 
hydrostatic equilibrium inside a dark matter halo potential; gas energy per unit mass similar 
to that of the dark matter, but modified by removal of low entropy gas and by feedback; 
appropriate outer boundary conditions; and pressure support from a non-thermal component. 
Some other processes necessary for detailed models will not be included because they are 
not in general required for obtaining global properties; though the results obtained here may 
need to be modified for those clusters having distinctly cooler cores (Allen & Fabian 1998). 
Finally, we will drop the limitation of a spherical NFW model and generalize to any case for 
which the dark matter potential is known. While several of the papers quoted above have 
allowed for some of these effects, none has included all in a fashion that can be adapted to 
an arbitrary gravitational potential. 

The next section reviews properties of the NFW model; §3 derives the properties of an 
initially parallel gaseous component and §4 derives these properties after the gas rearranges 
itself; §5 presents the resulting profiles and compares global properties with observed clusters. 
All of these sections assume spherical symmetry. In §6 we generalize the polytropic model 
in order to remove geometrical constraints, concluding our discussion in §7. 



2. The NFW Profile 

Navarro, Frenk & White (1997) have proposed a universal profile for dark matter halos, 
which we will first review here to establish nomenclature. Formally, the NFW profile extends 
to infinity and has logarithmically diverging mass; we will instead truncate the profile at the 
virial radius. In this section we first establish the properties of a distribution of matter with 
a truncated NFW profile. Assume that the density p depends on radius r as 

p(r) = I Kr+n) 2 r ^ r v , (1) 
[0 r > r v 

The virial radius r v is the radius within which the mean density is 200 times the critical 
density p c . The parameters C = r v jr\ and p\ define the model. It follows that the mass M 
is distributed inside r v as: 

M(r) = 47r p x r% g(x) , (2a) 

g(x) = ln(l+x)-— j-^ , (2b) 

where x = r fr\. Thus M tot = M(r v ) = An pir\ g(C) is the total mass. The rotation curve, 
V c 2 (r) = GM(r)/r = ATrGp 1 rfg(x)/x provides a useful label for the mass distribution; it has 
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a maximum at x 



c,max ~ 2.163, so the maximum circular velocity is 
V 2 _ A „n„Ji, 



c,max = ^Gp X r x G m{ 
G m ax — 9^Cc,max) / %c 



0.2162 



(3a) 
(3b) 



Note that a given M to t and V c>max can be used to completely define the matter distribution 
(as an alternative to C and p\). 



The gravitational potential from this mass distribution is 



where 



/(*) = 



G, 



ln(l+x) _ 1 
x 1+C* 
ln(l+C) _ _J_ 
C 1+C 



x<C 
£ x > C 

x 



The total gravitational energy is then given by 



w = - 



Mtot 



where 



H(C) = 



$(r)rfM = -G- m \XmaxM to tH{C) 
ln(l + C) C(l + C/2)" 



(i + c) (i + cy 



9(C) 



(4) 



(5) 



(6a) 



(6b) 



Assuming velocities are isotropic for the bulk of the matter distribution, the velocity 

iM 2 ) = ~Pt 



dispersion of the dark matter in 1-D, a 2 (r), obeys the equation -r-(pa 2 ) = — p4-&, which has 



the solution 

a 2 = ^S c (x)x(l+x) 2 , (7a) 



Gfria-j 



where 



S c (x) = S C (C) - ^ ^(1+^)3 ~ dx ' ( 7b ) 

with Sc(C) a positive constant, the value of which will be determined below. The corre- 
sponding pressure is simply P = pa 2 = AnGpfrfScix), so the average pressure is 



P = -j^-j [ Pdr = UnGplrlC- 3 [ S c {x)x 2 dx , (8a) 
and the surface pressure is 



P s = ATrGpirfSciC) . (8b) 
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A boundary surface pressure is required because of the jump in density and pressure at 
r = r v , and would be provided in any realistic physical model by the momentum flux from 
infalling matter at the boundary. One way to quantify the pressure term represented by 
Sc{C) is to estimate the momentum per unit area transported in by infalling matter. If the 
rate of mass accretion is M tot and the accreted mass (starting at rest from the turnaround 
radius) is moving at freefall velocity vj-j- = V^(r v ), then P s = M tot v f f / (4nr 2 ) . In the self- 
similar solution of Bertschinger (1985) M tot oc t 2 ^ 3 . This should be approximately correct, so 
we will take M tot = 2qM tot /(3t), with deviations allowed for with use of the correction factor 
q. Clusters of galaxies are by definition in overdense regions, so the surroundings will always 
appear to be close to the critical density p c , hence an appropriate time is t — (QnGp c )~ l l 2 ] 
thus we will adopt 

p s =4« v ^ps(- £ -) 1/2 > ® 

\PsPsJ 

where p s is the density at r v , and p is the mean density inside r v . Combining this with Eqn. 
(8b), it follows that 

o (n (2 2 9 3 (C) Pc \ 1/2 k/{C) (P\ 1/2 g\C) 

Sc{c) = \z q —Jj =V2q —{7j =q wc^ ■ (10) 

Alternatively, one may take the expression for the radial velocity dispersion in an NFW 
halo derived by Lokas & Mamon (2001) and evaluate it at the virial radius x — C. Lokas & 
Mamon (2001) solved the Jeans equation for an NFW density profile assuming a constant 
velocity anisotropy; for the simplest case of isotropic orbits, their Eqn. (14) yields 

Sc{C) = \ - - ^ - 2 (l+C) 2 ~ T+C + (I + 2^ ~ h ~ T+c) + ^) + (I]) 

|ln 2 (l + C) + 3Li 2 (-C) 1 J 

with the dilogarithm Li 2 (x) = J"°ln(l — y)d\mj. For the low concentration halos considered 
here, this gives pressures similar to Eqn. (10) with q m 4. We will use Eqn. (11) in the 
following, but have found simply using q = 4 gives very similar results. An intriguing feature 
of the NFW model is that the coarse-grained phase space density p(x)/a 3 (x) is a power law 
in radius, following x' 15 / 8 (Taylor & Navarro 2001; Williams et al. 2004). Using Eqn. (11) 
to set the surface pressure, the phase space density found by combining Eqns. (1) and (7a) 
does display in this power law behavior, with the appropriate slope. If we instead use (for 
example) q — 1, then the same power law still holds inside 0.5r v ; only nearer the surface is 
there a significant deviation. This insensitivity to the exact choice of surface pressure helps 
explain why the approximation q m 4 works well. 
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3. Initial Gas Energy and Pressure 

The goal of this paper is to populate with gas, in a physically plausible fashion, the 
potential well created by a dark matter halo with known mass, radius, concentration, and 
maximum circular velocity. Let the ratio of initial gas mass to total mass equal the cosmic 
average, i.e. the total gas mass M g ^ = Q b M tot /Q m . We begin by assuming that initially 
the two components have a parallel distribution — which is what would be expected in the 
absence of energy transport mechanisms. Thus the gas pressure is by hypothesis P g = 
(M g , t /M tot )P. 

Now we apply the Virial Theorem to the whole, with allowance for the non-negligible 
surface pressure: W + 2T — Anr^Ps = 0. Since the kinetic energy T is related to the mean 
pressure (in the absence of significant bulk motions) by 2T = Anr^P, we have, using Eqn. 
(8a), 

W = -Anr 3 v (P - P a ) = -3(a 2 )M tot (l - 8 S ) , (12) 

where 5 S = P s /P is the surface pressure in terms of the mean. When combined with Eqn. 
(6a), this becomes 

3(a 2 )(l-5 s ) = 0^X^(0) . (13) 

Note this is taken from the definition of (a 2 ) = j a 2 dM / M tot = (^nr^P) /M tot . Rewriting 
Eqn. (13) using Eqns. (6a) and (7a) gives the dimensionless form of the Virial Theorem: 

3(1 - S s ) [° S c (x)x 2 dx = H{C)g{C) . (14) 
Jo 

It follows that 

c _ S C (C) 

s S c (C) + C- 3 H(C)g(C) ' 1 ' 

Given S C (C), 5 S is known, and Eqn. (13) can be rewritten as (a 2 ) = K(C)V 2 max) where 

K{C) = \h(C)-0I- ) , (16) 

and the potential energy from Eqn. (6a) is now 

Wo = -3K(C)(1 - 5 s )M tot V 2 max . (17) 

Thus the kinetic and potential energies can be specified, once the NFW parameters for a 
dark matter halo have been specified. 



Earlier we assumed a monatomic gas component, initially distributed in the same man- 
ner as the dark matter. What is the total energy of the gas? Treating it as a tracer of 
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negligible mass, that is, assuming the gravitational potential is totally determined by the 
dark matter, E g = \M g>i (3a 2 } + 2 (M g>i /M tot ) W . Combining this with Eqn. (17) gives 

E 9 = -lM g ,V c 2 max K(C)(3 - A5 S ) . (18) 

Also, the gas surface pressure from Eqns. (8b) and (3a) is 

p _ M gji V 2 max S C (C) 
S ' 9as ~ G max g(C) ■ 1 ' 

Thus both the gas energy and surface pressure are now known in terms of the dark matter 
halo parameters. 



3.1. Allowance for Stellar Mass Dropout 

Star formation will change the amount of energy in the remaining gas, because that 
portion of the the gas which collapses and is removed will have a lower entropy and a shorter 
cooling time than is typical. In our idealized halo, this is the gas which would end up in 
the central region, so we will estimate the change in energy by removing a fraction of the 
core; this removed fraction (corresponding to the mass in stars) has the shortest cooling 
time and lowest entropy. Fukugita, Hogan & Peebles (1998) estimated that for clusters the 
stellar mass inside the virial radius is roughly 0.19/i ' 5 times the gas mass, and Sanderson & 
Ponman (2003) find a median stellar to gas ratio of 0.21; but lower values have been found 
by Balogh et al. (2001) and Lin, Mohr & Stanford (2003). We will adopt a stellar to gas 
mass ratio of 0.12 independent of cluster mass, in reasonable agreement with the results of 
the latter two papers, assuming a LCDM model with h=0.7. For simplicity we will assume 
enough gas is turned into stars for this ratio to hold generally. In other words, the ratio of 
the mass in stars to the final gas mass is f s = 0.12, with M g = M g ^j{\ + f s ) remaining in 
the gaseous state. This will be done by removing all gas inside a radius x s , found by solving 
9( x s) = g{C)f s / (1 + fs)- The initial energy in this remaining gas is found by integrating the 
previous expressions for the energy over only the mass M g outside of x s \ 

E 9 = + fs)M g V c 2 max [K(C)(3 - 45 s ) + K s (x s )/G max ] (20) 

with 



K s {x s ) 



9(C) 



Sc(x)x dx 



■2 , 2 f Xs f(x)xdx 



o 



(21) 



which, while unfortunately not as simple as before, can still be determined once the dark 
matter halo parameters are given. By removing this gas we are both lowering the gas mass 
and, more critically, increasing the energy per unit mass of the remaining gas (Voit & Bryan 
2001). 
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3.2. Other Changes to the Gas Energy 

The gas energy will change due to the the work done by any increase or decrease of the 
gas volume. To calculate this latter term, we will assume that the surface pressure does not 
change with radius, i.e. it is always given by Eqn. (19). Let Cf be the final outer radius of 
the gas distribution in units of r x ; then 

AEp = -f rl ^ - ° 3) Ps ^ s{c) = 4 (1+/ * )M ^— - c ^ ■ (22) 

In addition, we can expect there to be energy input to the gas from feedback processes. 
These are primarily of two kinds: wind and supernova shock energy deposited in the hot gas, 
and heating due to output from accreting massive black holes in the centers of the massive 
galaxies (Scannapieco & Oh 2004). Dynamical friction on galaxies moving through the gas 
at somewhat trans-sonic speeds may also be a source of energy input (Ostriker 1999; El-Zant, 
Kim & Kamionkowski 2004; Faltenbacher et al. 2004). To a good first approximation all of 
these effects are proportional to the gas mass (there is no evidence that the efficiency with 
which gas is transferred to stars varies strongly and systematically with V 2 max for moderately 
rich clusters) . Let e be a measure of the efficiency with which gas is heated by the condensed 
component; the energy input can then be written as 

AE f = ef s M g c 2 . (23) 

The value of e from energy output of supernovae can be estimated as the product of the 
fraction of mass turned into stars (/ s =0.12), the number of supernovae per solar mass ex- 
pected from a Salpeter initial mass function (0.007M" 1 ), and the energy input per supernova 
(10 51 erg); this gives e = 2.8 x 10 -6 or ef s = 3 x 10~ 7 . While assuming perfect efficiency 
in transferring this energy into the gas is unrealistic, we find that by itself this amount of 
energy has little impact on the gas profile, consistent with previous results (Balogh, Babul & 
Patton 1999; Valageas & Silk 1999; Bower et al. 2001; Shimizu et al. 2004). The energy input 
from AGN is more substantial. For a galaxy hosting a supermassive black hole, the ratio of 
black hole to stellar mass will be 0.0013 (Kormendy & Gebhardt 2000; Merrit & Ferrarese 
2001), so the ratio of black hole to total gas mass will be roughly 1.6 x 10~ 4 . Observational 
constraints give the efficiency with which energy is released from accreting black holes to be 
~ 0.10 (Yu & Tremaine 2002); adopting a conversion efficiency to mechanical energy of 30% 
(Inoue & Sasaki 2001) leads to an efficiency ef s = 4.7 x 10" 6 or e = 3.9 x 10~ 5 . We will 
adopt this last number in the rest of the paper. This is equivalent to an energy input of 2 
keV per particle or 4 keV per baryon. Another way of looking at this is to divide ef s M g c 2 by 
a Hubble time to compute a typical luminosity; for a halo with a total mass of a few times 
1O 14 M , this is of order 10 44 erg/s. In other words, the energy input from black holes is the 
same magnitude as the observed energy radiated in X-rays. 
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4. Polytropic Rearrangement 

Now assume that the gas rearranges itself (changing its density profile) through un- 
specified processes into a polytropic distribution with polytropic index n or adiabatic index 
T = 1 + n~ l . The commonly addressed cases are (r, n) = (1, oo) for isothermality, and 
(5/3,3/2) or (4/3,3) for non-relativistic or relativistic isentropic fluids respectively We know 
that there is turbulence (Kim & Narayan 2003a) induced by merger shocks (Bryan & Nor- 
man 1998) and galaxy wakes (Stevens, Acreman & Ponman 1999; Sakelliou 2000). Other 
processes driving the rearrangement could include, for example, conduction (Kim & Narayan 
2003a; Dolag et al. 2004) or wave transport of energy (via gravity, Alfven, or magnetosonic 
waves, e.g. Cen 2005). This rearrangement means that the outer gas radius could be larger 
or smaller than r v . 

For central pressure Po and density po, a polytropic distribution requires that the gas 
pressure P' and density p' after rearrangement are related by 

^' = W/Po) (1+1/n) = W/Po)' , (24) 

and the central isothermal gas sound speed is defined as V^ Q = Pq/ p' . Note that T is not 
in fact the actual ratio of specific heats; we require only that the gas has arranged itself in 
polytropic fashion, as in Eqn. (24). Fig. 1, kindly provided by Greg Bryan, shows results 
from a high resolution adiabatic AMR simulation of a massive (~ l^h^M®) cluster. The 
pressure-density relation in this calculation closely fits that of a V = 1.15 polytrope, and 
lends support for the proposal that turbulent mixing (the only one of the processes listed 
above that was included in this computation) can lead to a fairly tight polytropic relation. 
SPH simulations by Lewis et al. (2000) resulted in a pressure-density relation well described 
by a polytropic equation of state with T « 1.2; similar results are reported by Ascasibar 
et al. (2003) and Borgani et al. (2004). This result holds for both adiabatic and radiative 
simulations (but see Kay et al. 2004), and agrees well with the effective T derived from 
observed clusters by Finoguenov, Reiprich, & Bohringer (2001). Solanes et al. (2005) find 
that T = 1.2 offers the best consistency with the assumption that the specific energy of the 
hot gas equals that of the dark matter. Interestingly, the purely adiabatic, spherical, and 
self-similar collapse solution of Bertschinger (1985) was also polytropic with T w 1.17. 

An additional contribution to the pressure may come from a relativistic component; 
such a component could be created for example at shock fronts, converting part of the gas 
energy into cosmic rays. While the relativistic portion of the gas will contain a negligible 
fraction of the mass, it may contribute significantly to the total gas energy and pressure. 
We allow for the fact that, in addition to the gas pressure P' there may be a nonthermal 
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component having pressure P re \ = 5 re iP' with total pressure 

P tot = P' + P rel = P'(l + 5 rel ). 

It is not obvious if S re i is maximum in the center, where there may be injection of a relativistic 
fluid by an AGN, or in the outer parts of the cluster, where there may be injection of 
relativistic particle energy in boundary shocks. Thus, for simplicity we will take 5 re i = 
constant. 

Given these relations, the equation of equilibrium for a spherically symmetric distribu- 
tion, dPtot/dr = —pd&/dr, becomes 



d d f d\ il+1/n) 



GM tot (r) 



Thus 



ft 



A) V&(l + n)(l + 6 rel ) (l + n)(l + 8 rel ) 

This last is from Eqn. (4) with 



= AltGpirl = Vc,max r _i 



V 2 V 2 



(25a) 



and 



Thus, 



where 



x < C 

J{X) = ^ + x >c . (25b) 



P'{r) = p' 



[l + n)(l + 5 rel ) 



PT , (26a) 



^- (irf^M (26b) 

is the familiar polytropic variable defined by Chandrasekhar (1967). Eqn. (26a) was first 
derived by Makino, Sasaki & Suto (1998) for the isothermal case, and Wu, Fabian & Nulsen 
(2000) more generally. The final gas radius can be smaller or larger than the initial value. 
Denoting this radius in units of r\ as Cf, the gas mass can be written M g = 4nrfp' L, where 

L = L(n,P,C,C f )= 9 n x 2 dx . (27) 

Jo 
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The thermal component will contribute a factor of | J P'd 3 x to the kinetic energy, and 
the relativistic component 3 J 5 re iP'd 3 x, so the rearranged total gas energy is E' — |(1 + 
25 rei ) J™ 9 V s 2 dM + J^ 9 <$> tot dM. Defining two more integrals 



we now have 



as the final energy. 



h = I 2 (n,(3,C,C f )= / f(x)6 n x 2 dx 

Jo 

h = h(n,/3,C,C f ) = [ 1 6 1+n x 2 dx 



E' g = -M g V s 2 



(28a) 
(28b) 

(29) 



5. Constraints on the final temperature 

Suppose we have a dark matter halo for which the relevant properties — M tot , r v , C, V c>maxi 
etc. — are known. From the previous section, the final distribution of the gas can be de- 
termined as a function of the two unknowns (5 and C/; thus to specify the final gas tem- 
perature and density distribution it remains only to constrain these two parameters. The 
first constraint is from conservation of energy: the final gas energy will equal the initial 
energy plus changes to due star formation, expansion or contraction, and feedback; i.e. 
E g + AE P + AE f = E' g . Combining this with Eqns. (20), (22), (23), and (29) yields 

|(1 + fs) [G max K(C) (3 - 45 s ) + K s (x s )} - G max ef sV ^- 

* v c,max I On \ 

(keeping in mind that L, J 2 , and I3 are functions of j3 and C/). A second constraint comes 
from the fact that the surface pressure of the gas must match the exterior pressure, which 
we have fixed at P s , gas (C). This gives 



:i + f s )^-PL = (i+5 rel ) 



mcf) 



(l + n)(l + 5 rel ) 



l+n 



(31) 



Thus, given r v , C, and V c ^ max for the dark matter halo, and appropriate choices for T, 
f s , e, and 5 re i, Eqns. (30) and (31) can be solved for (5 and C/, and the central temperature 



(32) 
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is known, as is the density parameter p' Q , so the gas distribution is fully specified. The 
expected X-ray luminosity Lx can then be calculated. Following Balogh, Babul & Patton 
(1999), we will include both Bremsstrahlung and recombination, which becomes important 
for temperatures below 4 keV, by using the cooling function 

A(T) = 2.1 x lO-^T 1 / 2 [1 + (1.3 x 10 6 /T) 3 / 2 ] cm^rgs" 1 . (33) 

5.1. Simulated Halo Catalogue 

The plausibility of this procedure can be evaluated by trying it out on a population 
of many dark matter halos and comparing the results to observed clusters. The halos we 
use here come from an N-body simulation designed to be in concordance with observational 
constraints. The simulation is of a periodic cube 1500/?, _1 Mpc on a side containing N= 
1260 3 = 2 x 10 9 particles. The cosmology was chosen to be a standard LCDM power 
law model with the following parameters: baryon density = 0.047; Cold Dark Matter 
density ^cdm — 0.223 (hence total matter density Q m = 0.27); cosmological constant Q\ = 
0.73 (thus spatially flat); Hubble constant given by h — iJ /(100km s _1 Mpc _1 ) = 0.70 
(hence Vl^ti 1 = 0.02303); primordial scalar spectral index n s = 0.96; and linear matter 
power spectrum amplitude = 0.84. These values are consistent within one standard 
deviation to those derived either from WMAP data or from WMAP combined with smaller 
angular scale CMB experiments and galaxy data (Spergel et al. 2003). The initial conditions 
were generated using the publicly available code GRAFIC2 (Bertschinger 2001) to compute 
initial particle velocities and displacements from a regular grid. Since the memory required 
to hold a 1260 3 grid is 8 gigabytes, it was necessary to modify the single level portion of 
this program by adding message-passing commands in order to distribute the mesh among 
several processors. 

The simulation was carried out with the TPM (Tree-Particle-Mesh) code (Bode & Os- 
triker 2003), using 420 processors on the Terascale Computing System at the Pittsburgh 
Supercomputing Center; it took not quite five days of actual running time. The box size and 
particle number determine the particle mass of 1.26 x l0 u h~ 1 M Q . The cubic spline softening 
length was set to 17/i _1 kpc. A standard friends-of-friends (FOF) halo finding routine was 
run on the redshift z = box, using a linking length b = 0.2 times the mean interparticle 
separation (Lacey & Cole 1994); this yielded 575,125 halos with both a FOF mass above 
2 x lO 13 /i _1 M and a virial mass above 1.75 x 10 13 /i^ 1 M o . The PM mesh used in TPM con- 
tained 1260 3 cells, and at redshift zero all PM cells with an overdensity above 39 were being 
followed at full resolution, so these objects had the full force resolution of TPM. For the 
range of parameters used here, clusters with kT > 2keV contained more than 200 particles 
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within r v . 

For each halo, the position of the most bound particle is taken to be the cluster center. 
Then M tot and r v are measured, as are V c ^ max and the radius of maximum circular velocity 
r c ,max] this latter gives the concentration C = 2.163r v /r c ^ max . This defines the equivalent 
NFW model halo, i.e. that NFW model closest to the computed dark matter halo. With 
this information, the procedure outlined above can be carried out on each halo to compute 
the gas density and temperature. 

5.2. Resulting Profiles 

Given this set of halos, it remains to specify T, f s , e, and 5 re i. Let us first consider 
the appropriate T. Fig. 2 shows the projected temperature profile for different values of 
T; this profile was computed by integrating the emission- weighted temperature along the 
line of sight. To normalize the curves, the mean temperature (T) was calculated by evenly 
weighting all radii inside r v /2; this was done to correspond with the method of De Grandi 
& Molendi (2002), who measured the mean profile for clusters with and without cooling 
flows — shown in the Figure as filled and open circles, respectively. Examination of Fig. 2 
shows that r=l. 2-1.4, corresponding to polytropes with index n— 2.5-5, provides adequate 
fits to the outer parts of the clusters, within which resides most of the gaseous mass (see also 
the discussion in Solanes et al. 2005). Ascasibar et al. (2003) have shown that a V — 1.18 
model is a good fit to the average temperature profile measured by Markevitch et al. (1998); 
this latter measurement has been confirmed by De Grandi & Molendi (2002), Piffaretti et 
al. (2005), and Vikhlinin, et al. (2005). As discussed above (§4), T = 1.2 is also a good fit 
to hydrodynamical simulations (Lewis et al. 2000; Loken et al. 2002; Ascasibar et al. 2003; 
Borgani et al. 2004; Kay et al. 2004). The lack of an isothermal core will not lead to a 
serious overestimation of luminosity or emission-weighted temperature because the volume 
of this central region is small. This was tested by taking the T=1.2 profile and imposing an 
isothermal core, matching the density and pressure at 0.2?^; the resulting changes in emission- 
weighted temperature and luminosity were generally less than 10%. However, neglecting 
cooling will reduce the scatter in the M 500 — T and L x — T relations (McCarthy et al. 2004). 

The effect of the polytropic rearrangement on the radial profile of the gas can be seen 
in Fig. 3. The example halo used here has physical parameters M tot = 4 x 1O 14 /i^ 1 M , 
r v = 1.2/i _1 Mpc, C = 4, and V^ max = 1200km/s. The top two panels show the temperature 
(relative to T ew , the mean emission- weighted temperature inside a radius -R500 containing an 
density of 500p c ) and density (relative to 200p c ) as a function of radius. It is instructive 
to compare to the original NFW distribution, shown as a dot-dashed line; in this case 
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the central temperature goes to zero, as the density profile has a cusp. The polytropic 
rearrangement (taking T = 1.2 and f s — e — 5 re i = 0, shown as a dotted line) increases the 
central temperature while decreasing the density, removing the cusp (with a correspondingly 
dramatic lowering of the X-ray luminosity, as we shall see). This is seen more clearly in the 
third panel, which shows the ratio of gas to dark matter mass interior to a given radius, in 
terms of the cosmic average: inside the dark matter core radius, the gas fraction declines 
sharply. These temperature and density profiles result in the "entropy" profile shown in the 
final panel of Fig. 3, taking the definition of entropy to be Tp~ 2 / 3 . The polytropic profile 
has a slope close to r 11 near the virial radius, and is shallower nearer the cluster core; 
this behavior has in fact been observed in a wide range of clusters (Ponman, Sanderson 
& Finoguenov 2003; Pratt & Arnaud 2005; Piffaretti et al. 2005). This behavior has been 
derived before in analytic models assuming the gas is shock heated (Tozzi & Norman 2001), 
and is also seen in hydrodynamic simulations (Lewis et al. 2000; Borgani et al. 2004; Kay et 
al. 2004). 

This change in profile has a strong impact on other observable cluster properties, as is 
shown in Fig. 4. In these plots the temperature is taken to be the mean emission-weighted T 
inside -R50O) T ew . In clusters with more complicated structure this measure may not coincide 
well with the spectroscopically measured T, as pointed out by Mazzotta et al. (2004), who 
provide an alternative measure. However, for the simplified models here, the difference 
between emission weighting and the Mazzotta et al. (2004) spectroscopic-like measure is only 
a few percent at most. The lines show the median value as a function of T ew , found using the 
dark halo catalogue described in §5.1. The first impact of the polytropic rearrangement is 
to increase the observed temperatures. This is clearly demonstrated in the left-hand panel, 
which gives the mass-temperature relation. The points are from Reiprich & Bohringer (2002), 
as adjusted by McCarthy et al. (2004); here the mass is M 50 o, the mass inside a sphere 
containing mean density 500p c . The polytropic model distribution (dotted line) resembles 
that assuming an NFW profile, only shifted to higher T; note assuming an NFW gas profile 
leads to significant disagreement with the observed relation, while switching to a polytropic 
model provides much superior agreement. This is also true in the right-hand panel, which 
shows the bolometric X-ray luminosity as a function of T; the data points are the subset 
of the ASCA cluster catalogue (Horner 2001) described in McCarthy et al. (2004). The 
polytropic model, without the central cusp, yields a lower luminosity than the NFW profile. 
The slopes of both the L x — T and M 50 o — T relations retain the same self-similar values, 
however. 

The next physical input is the fraction of gas which collapses into stars. As discussed 
above, this is roughly one eighth the the gas mass inside the virial radius, or f s = 0.12. Since 
the stars in clusters are old, this fraction will hold for all moderately low redshifts. As shown 
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as short-dashed lines in Fig. 3, assuming f s = 0.12 for a typical cluster (keeping T = 1.2) 
increases the temperature slightly and reduces the gas density. Since the temperature change 
is not large, this has little effect on the M 50 o — T relation. However, gas removal for star 
formation, which increases the mean energy per particle for the remaining gas, leads to 
lower densities and so has a significant impact on the L x — T relation. For the most massive 
(hottest) clusters, the predicted luminosity is in fact close to that observed; however, the self- 
similar slope of the relation is still preserved, so for less massive clusters L x is overestimated. 

The next required physical input is the amount of energy from feedback coming from 
supernovae and active galactic nuclei, discussed in §3.2. The results of including feedback 
of e = 3.9 x 10~ 5 in the V = 1.2, f s = 0.12 model are shown as long-dashed lines in Figs. 3 
and 4. As one would expect, the radial profile has a higher temperature and lower density. 
However, the effect of feedback differs from those considered previously, because the resulting 
relations are no longer self-similar. For massive clusters with V c>max > lOOOkm/s or kT > 
lOkeV, feedback is of little importance because the added energy is small compared to the 
gravitational energy, but for smaller masses it can have a significant impact. One can see a 
steeper slope in the M500 — T relation, but the most significant effect is on the luminosity, 
which in shape now more closely resembles the observed distribution. 

The remaining physical effect left to include is nonthermal pressure. We will take 
8 re i = 0.1; the nonthermal sources of pressure may in fact contribute a few tens of percent 
of the total (Miniati 2004). The results can be seen by comparing the solid (5 re i = 0.1) and 
long-dashed (5 re i = 0) lines in Figs. 3 and 4. With this additional support, less kinetic energy 
is required at a given pressure. Thus, while the density profile is little changed, the resulting 
gas distribution is somewhat cooler, and the emission weighted temperature is lower at a 
fixed M 500 or L x . 

The departure from self-similar scaling is shown further in Fig. 5, which displays the 
radial profiles of temperature, gas density, gas fraction, and entropy for clusters of mass 
M tot = 10 15 ,5 x 10 14 ,2.5 x 10 14 ,1.25 x 10 14 , and 6.25 x IO^/i^Mq. Star formation and 
feedback were included with f s = 0.12 and e = 3.9 x 10~ 5 , but not a relativistic component. 
For the least massive cluster we took C — 4, V CjTnax = 700km/s, and r v = 650/i _1 kpc, scaling 
for the others as C oc M~°- 13 , V c>max oc M 1//3 , and r v oc M 1 / 3 ; this is in reasonable agreement 
with our N-body cluster catalog. With these parameters, the emission weighted temperatures 
inside i?soo are kT ew — 10.1, 6.7, 4.5, 3.1, and 2.2 keV, respectively. For decreasing mass, the 
temperature and density profiles become increasingly shallow, leading to a faster decrease 
in X-ray luminosity. The ejection of gas following feedback energy injection leads to a gas 
fraction (relative to the universal value) less than unity at the virial radius; with the full 
halo catalog and these parameters we find for halos in the range 1 — 2 x lO 14 /i _1 M , the gas 
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fraction at the virial radius is 0.72±0.09 (one standard deviation). This result is in agreement 
with simulations including both heating and cooling: Muanwong et al. (2002) and Kravtsov, 
Nagai & Vikhlinin (2005) find that for halos with M ~ 1O 14 /i _1 M the hot gas fraction 
inside a radius enclosing overdensity ~ 100p c is in the range 0.6-0.7, while Ettori et al. 
(2004) find slightly higher values of 0.7-0.8. Observational estimates give similar values with 
a higher scatter (Evrard 1997; Mohr, Mathiesen & Evrard 1999; Sanderson & Ponman 2003). 
Both these observational and the theoretical studies suggest that more massive clusters have 
higher hot gas fractions, behavior which is reproduced here. 

The bottom panel in Fig. 5 displays the entropy T/nj^ 3 , with n e the electron density 
in units of cm -3 (and T in keV), which can be compared directly with the observations 
listed above. The entropy profiles have a slope close to the observed value of r 11 at the 
virial radius, but this slope quickly becomes shallower for smaller radii. However, we have 
not taken into account cooling, which would be more important near the center, steepening 
the inner entropy profile (McCarthy et al. 2004). It has been observed that the entropy 
in clusters scales as T a65 (Ponman, Sanderson & Finoguenov 2003; Pratt & Arnaud 2005; 
Piffaretti et al. 2005), rather than linearly in temperature as one would expect for self-similar 
scaling. Both of these scalings are shown in Fig. 6; the upper panel gives T jnj z jT ew and 
the lower panel T jnj 3 /T^ 65 . It is clear from this Figure that the observed T 65 scaling is 
followed quite closely. However, at radii near r v the observational picture is unclear; with 
a sample of 14 nearby clusters, Neumann (2005) found that the outer regions followed self- 
similar scaling and may be affected by the accretion of cooler material. Also, for poorer 
systems than considered in this paper (kT < 2keV), we find that this scaling breaks down; 
such groups may have different entropy profiles than are seen in richer systems (Mahdavi et 
al. 2005). 



6. Generalization to an Arbitrary Dark Matter Potential 

The virtue of the method presented in this paper is not just in the equations presented 
in the preceding sections, whose purpose was to establish physical and mathematical prin- 
ciples and assess the plausibility of the results, but also in its ability to efficiently model 
complex asymmetrical systems containing substantial substructure. In this section we relax 
the previous assumption of spherical symmetry and apply the same method to more complex 
DM potentials. Suppose that a cluster potential is known from an accurate DM integration; 
the cluster will likely be aspherical and contain significant substructure. Then, if we are 
satisfied by the ability of an equilibrium polytrope to model the gas in a cluster of galaxies, 
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the integration of the equation of equilibrium VPm = — pV0 gives 

(1 + 5 rel ) F P ^ = -0(r) + (1 + Wp^r- + 0o • (34) 

r - 1 p(f) r - 1 p 

The last two terms comprise a constant of integration; here O is the potential minimum 
located at position r = fo, and the pressure and density at this point are designated by P 
and pq. Then, making the definition 

d{r) = 1 + h T ~ 1 ^ M> - , (35) 

(1 + d re i)T P 

the pressure and density are simply 

P(f) = P 9(r)^i (36a) 
p(f) = p 6(r)^ (36b) 



where 9{r) is essentially the same polytropic variable defined by Chandrasekhar (1967). 
Thus for an equilibrium polytropic gas residing in a known potential (f>(r), the determination 
of the structure is reduced to the determination of the two numbers Po and p . Adopting 
the approach taken in the previous sections, these constants can be determined by satisfying 
two equations of constraint on the final energy and the surface pressure. 

We carried out this procedure on the same N-body simulation used previously in the 
following manner. A set of particles identified as a cluster is placed in a nonperiodic 3-D grid. 
The grid cell size I is set to four times the N-body particle spline softening length, as scales 
smaller than this can be affected by numerical resolution issues; increasing or decreasing I by 
a factor of two had little impact on the results. The dark matter density in each cell k, pk, is 
found from the particle positions using cloud-in-cell (CIC), and the gravitational potential <pk 
on the mesh is calculated from the density using a nonperiodic FFT (Hockney & Eastwood 
1981). The position of the cell with the minimum potential MIN(<f)k) = 4>o is taken to be 
the center of the cluster, rj. The cluster velocity is estimated as the mean velocity of the 125 
particles closest to r^; this mean is subtracted from all the particle velocities. Then the DM 
kinetic energy per unit volume tk = pv 2 is found in each cell: as with the mass, the kinetic 
energy of each particle is distributed among 8 cells using CIC. The virial radius r v and DM 
mass M to t are found from the density distribution. The N c i cells inside r v are identified, as 
are the JV& cells in a buffer region of width r^/ surrounding the cluster with centers in the 
range r v < r < r v + r bu f. The buffer width was set at 9 cells, or r bu j = 153/i _1 kpc for the 
simulation used here. The gas surface pressure is taken to be the mean value (assuming 
velocities are isotropic) inside this buffer region: 

P - = V£^T** (37) 
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where the sum is over all cells in the buffer region. Assuming gas traces the DM, the gas mass 
inside r v is originally M 9ji = ^-M tot . As before, the portion of this gas which is turned into 
stars is f s M 9ti /(l + f s ). To decide which portion of the gas becomes stars, cells are ranked 
by binding energy Pk4>k + \tk\ starting with the most bound cell, the initial pk and i& are set 
to zero for each cell in turn until the gas mass removed, ^p^l 3 , totals to f s M gji /(l + f s ). 
The original gas mass inside r v is then 

the sum being over all cells inside r v , and the initial gas energy is 

Ncl n 

^ = E?r{^ + ^}^ • ( 39 ) 

fe=l 

As before, this energy can be supplemented by feedback energy AEf = ef s M g c 2 . 

As in the 1-D case, this gas is assumed to rearrange itself into a polytropic distribution 
with T = 1.2. It only remains to specify Pq and po, which are fixed by the final energy and 
surface pressure. For a given initial choice of (P ,p ), the final gas density and pressure 
can be found after calculating 6 k for each cell from Eqn. (35). As before, the initial energy 
may be changed by the inflow or outflow of gas. The final radius of the gas, 77, is found by 
moving outwards from the cluster center until mass M g is enclosed, i.e. 

Nj 

M g = X>o0r^ 3 , (40) 
fc=i 

where the sum is over the Nf cells inside 77. Similarly to Eqn. (22), we will assume that the 
surface pressure does not change with radius, so the change in energy due to expansion or 
contraction is proportional to the change in volume. This means AEp = (47r/3)(r^ — r 3 f)P s , 
with P s given by Eqn. (37). Now we have all the information required for the first constraint 
on (P ,p ), namely the conservation of energy, where the final gas energy is 

Ef = J2 {poC^o + )l 3 = E 9 + AEf + AEp . (41) 



k=i 



The second constraint is the mean pressure in the N^j buffer cells between 77 and 77 + r\, u f\ 
this is assumed to match the original value: 

N b :}J2 P ^=Ps ■ (42) 
k=l 
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Thus after an initial estimate for (P , p ), it is now possible to iterate to a solution satisfying 
Eqns. (41) and (42). This solution provides the full three dimensional pressure and density 
of the gas with allowance for substructure, triaxiality, etc. 

An example of the resulting gas distribution for one halo taken from the catalog is given 
in Fig. 7. The simulation particle positions are shown in the upper left-hand panel; there 
are several large substructures in the process of merging with the main objects. The volume 
shown is a cube 6.4/i _1 Mpc on a side. The upper right-hand panel shows the projected gas 
surface density obtained by the method just described. With the gas density and temper- 
ature, maps can be made for the X-ray emission and SZE, as shown in the lower panels. 
The scale is linear, with black a factor of 100 below white. The results of this procedure 
employed on the entire set of halos used before are displayed in Fig. 8. A stellar fraction 
f s = 0.12 and feedback e = 3.9 x 10 -5 were included. For each plot, the median value is 
shown as a solid line and the shaded region encloses 68% of the halos at that temperature. 
Also shown as dashed lines are the results of the method of §5 based on the NFW profile; the 
scatter seen using this latter method is somewhat smaller than that shown for the full 3-D 
method, as the latter realistically includes substructure and triaxiality. We are neglecting 
cooling, which would increase the scatter further, and tend to increase the luminosity at 
a given temperature (McCarthy et al. 2004). Use of the NFW approximation appears to 
have little effect on either the M 500 — T or the L x — T relation, relative to using the full 
particle distribution. Examination of individual clusters shows that the spherically averaged 
gas profiles resulting from the N-body potentials are slightly shallower in both temperature 
and density, although the amount of gas inside R 500 is the same for both methods; the larger 
clumping factor when using the true density, with triaxiality and substructure, increases the 
luminosity enough to compensate for this. 



7. Conclusion 

The NFW model has provided a useful description of the distribution of matter inside 
collisionless DM halos, such as those hosting X-ray clusters. This success inspires hope that a 
similarly concise description can be found for the hot gaseous component. Given a population 
of dark matter halos from an N-body simulation (or from some semi-analytic model such as 
extended Press-Schechter), can we deduce the global properties of the baryonic component 
inside each halo? In this paper we have worked towards providing a prescription which is 
simple enough to apply broadly while remaining physically well motivated. 

In the model presented here, the gas is assumed to initially have energy per unit mass 
equivalent to that of the dark matter; this energy can be modified by removal of low entropy 
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gas (to form stars), addition of feedback energy expected from supernovae and accreting 
black holes, and mechanical work done as the gas expands or contracts. The gas is assumed 
to redistribute itself into a polytropic distribution in hydrostatic equilibrium with the DM 
potential; given the constraints on the total energy and the surface pressure at the virial 
radius, the gas distribution is entirely specified. 

We applied two variations of this method to a catalogue of cluster-sized dark matter 
halos drawn from a large cosmo logical N-body simulation. In the first variant, the mass, virial 
radius, and concentration of each halo was measured, and the mass profile was assumed to 
follow a spherically symmetric NFW profile. To determine the gas distribution in this case 
means solving Eqns. (30) and (31). We then allowed for complex, nonspherical profiles and 
substructure by using the full set of particle positions and velocities in each N-body halo 
to determine the potential and kinetic energy. These two methods give similar results, but 
assuming a spherical NFW profile gives slightly lower temperatures on average and gives 
significantly less scatter in the mass-temperature and X-ray luminosity-temperature than is 
observed. 

Simply assuming the gas follows the dark matter leads to too low temperatures and 
too high central densities, since the DM profile has a cusp. The polytropic rearrangement 
increases the central temperature while decreasing the density, removing the cusp. Remov- 
ing low entropy gas for star formation further increases temperatures and reduces density. 
However, neither of these processes changes the self-similar nature of the model. Including 
energy from feedback does change this, because in massive clusters the energy input will be 
small compared to the total binding energy, while for smaller masses it can have more of 
an impact. We also implemented a simple approximation for including nonthermal pressure 
support; including a relativistic component in this way leads to somewhat lower temperatures 
and slightly higher densities. 

Essentially two dimensionless numbers are required to prescribe the state of the gas in a 
given DM potential: the fraction of gas mass transformed to a condensed (primarily stellar) 
form (determined by observations to have the value f s ps 0.12); and the feedback from the 
condensed component, for which a plausible estimate for the energy output from supernovae 
and black holes that is trapped in the cluster gas is ef s ^ 0.12 • 3.9 x 10 -5 = 4.7 x 10~ 6 . 

The utility of fully understanding the properties of the intergalactic medium in clusters 
can be seen in Fig. 9, which shows the cumulative temperature function; the data points 
are from Ikebe et al. (2002). The lines (the line types are the same as in Fig. 3 and Fig. 4) 
demonstrate how, as different processes are included, the resulting temperature function can 
change quite dramatically. (Note that the curves are the number density at z — 0, whereas 
at the highest T objects are sufficiently rare such that the data points actually reflect the 
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cluster density at z > 0, which is lower). There is an apparent conflict in that no single 
model seems to satisfy all the observations simultaneously but this may not be a serious 
problem, and may in fact reflect the strength of the models. In order to reproduce the 
observed M 50 o — T and L x — T relations requires a high T and low density, in other words 
a significant amount of star formation and feedback. But this seems to predict too high a 
number density of clusters for a given T, as is seen in the plot of n{> kT). However, rather 
than a problem with the gas model, this may simply be due to our choice of cosmological 
parameters Vt m and a 8 when generating the cluster catalogue. The mass function from our 
N-body simulation is also too high when compared to that observed by SDSS (see Fig. 2 of 
Younger, Bahcall & Bode 2005). Thus an accurate mass-temperature relation would also 
lead to an overestimate in the temperature function, so the failure to reproduce the observed 
n(> kT) relation reflects a failure of the cosmological model; lowering og and/or Q m would 
alleviate this problem without significantly altering the predicted M 50 o — T and L x — T 
relations. A 10% reduction in a§ would reduce the number of clusters with kT > 4 by 
roughly a factor of two. This points out the usefulness of the cluster number density as a 
probe of cosmological parameters, but also the necessity of including all the relevant physics 
accurately. 

A variety of telescopic surveys in many wavelength bands will soon greatly multiply 
the number of galaxy clusters catalogued, particularly at higher redshifts. Unlocking the 
power of these new observational datasets as cosmological probes will require sophisticated 
theoretical predictions. In the future we plan to apply the methods developed here to explore 
the properties of clusters at higher redshifts and make detailed predictions for X-ray and 
SZE surveys in many different cosmological models. 
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Fig. 1. — Pressure as a function of density, from an adiabatic simulation of a massive cluster 
by Greg Bryan. The line has a logarithmic slope of T=1.15. 



-27- 




0.2 0.4 0.6 0.8 1 

R/r v 

Fig. 2. — Projected temperature profile for T =1, 1.1, 1.2, 1.3, 1.4, and 1.67 (solid lines, 
from top to bottom at outer radius). < T > is calculated by evenly weighting all radii 
inside r v /2. The points are from De Grandi & Molendi (2002); filled and open circles are 
clusters with and without cooling flows. The parameters used to calculate the profile are 
C = A,V Cjmax = 1200 km/s, M tot = 4 x lO 14 /i~ 1 M , and r v = 1.2/i -1 Mpc. Aside from the 
innermost region, the model with a T of 1.2-1.4 satisfactorily represents the T falloff in the 
outer regions, which contain most of the gas mass. 



-28- 




^^"'"f '"l I I I I I I I I I I I I I I I |_^ 

0.2 0.4 0.6 0.8 1 



Fig. 3. — Radial distributions of the temperature, density, mass fraction, and entropy of 
the gas for different model assumptions. Dot-dashed: gas follows NFW profile. Dotted: 
hydrostatic equilibrium with polytropic index T = 1.2. Short-dashed: also removing 12% 
of the gas to form stars. Long-dashed: also including energy input from feedback. Solid: 
also with a relativistic component. In the bottom panel, the entropy is scaled by a factor of 
T ew /(200p c ) 2 / 3 . 
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Fig. 4. — Results of the model applied to a population of DM halos. Line types as in Fig. 3; 
lines show the median value as a function of emission- weighted T inside overdensity 500. 
Left: mass (measured at overdensity 500). Right: bolometric X-ray luminosity. The data 
points in these two plots are from McCarthy et al. (2004). 



-30- 



i i i r 




0.1 1 

r/r v 

Fig. 5.— Radial distributions for masses M tot = 10 15 ,5 x 10 14 ,2.5 x 10 14 , 1.25 x 10 14 , and 
6.25 x 10 13 /i _1 M o (line types from top to bottom in uppermost panel). Shown, from top, are 
temperature, density relative to 200 times critical, the gas mass fraction interior to r, and 
the entropy T jnj z (here n e is in units of cm -3 ). Parameters used are T = 1.2, f s = 0.12, 
e = 3.9 x lO" 5 , and 5 rel = 0. 
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Fig. 6. — Radial distribution of entropy T jnj^ for cluster masses M to t = 10 15 , 5 x 10 14 , 2.5 x 
10 14 , 1.25 x 10 14 , and 6.25 x 1O 13 /i _1 M (line types as in Fig. 5) scaled by T ew (top panel) 
and by T^ 65 (bottom panel). 



-32 - 




Fig. 7. — Adding gas to a dark matter halo (drawn from a larger simulation) containing 
significant substructure. Upper left: DM particle distribution. Upper right: projected gas 
surface density. Lower left: X-ray luminosity. Lower right: integral of the pressure along 
the line of sight, proportional to the SZE decrement. The color scale is normalized to the 
maximum value; the volume shown is 6.4/i -1 Mpc on a side. 
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Fig. 8. — Results from using the directly computed three-dimensional potential of the simu- 
lated halos. The solid line in each panel is the median value, and the shaded region encloses 
68% of the clusters. The corresponding NFW approximation is shown as a dashed line. 
Parameters used are T = 1.2, f s = 0.12, e = 3.9 x 1CT 5 , and 5 re i = 0. Points are the same as 
in Fig. 4. 




Fig. 9. — The cumulative temperature distribution function, varying parameters as in Figs. 3 
and 4. (using the same line types); data points are from Ikebe et al. (2002). This function 
will be very sensitive to the choice of cosmological model, particularly Q m and cr 8 . 



